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Abstract 

The semiclassical SU(3) Skyrme model is traditionally considered as describing a rigid 
quantum rotator with the profile function being fixed by the classical solution of the cor- 
responding SU(2) Skyrme model. In contrast, we go beyond the classical profile function 
_ by quantizing the SU(3) Skyrme model canonically. The quantization of the model is per- 

formed in terms of the collective coordinate formalism and leads to the establishment of 
purely quantum corrections of the model. These new corrections are of fundamental im- 
portance. They are crucial in obtaining stable quantum solitons of the quantum SU(3) 
Skyrme model, thus making the model self-consistent and not dependent on the classical 
solution of the SU(2) case. We show that such a treatment of the model leads to a family 
of stable quantum solitons that describe the baryon octet and decuplet and reproduce the 
experimental values of their masses. 
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1 Introduction 



The Skyrme model is a nonlinear field theory having localized solutions, the so-called skyrmions, 
that are of finite energy and are characterized by a topological charge. It is an effective theory 
of low-energy QCD in the limit of a large number of colours, thus it describes baryons in a 
weakly coupled phase as was initially argued by [HE1E]. Indeed, the semiclassical quantization 
of the model has proven to be successful in describing the phenomenological properties of the 
baryons in the low-energy region. 

The SU(2) Skyrme model was originally defined to describe a unitary field U(x, t) in a 
fundamental representation of the SU(2) group with a natural boundary condition U — > 1 at 
the spatial infinity, |x| —> oo. This implies that the unitary field represents a topological map 
S 3 — > S 3 with an integer-valued winding number classifying the solitonic sectors of the model. 
This topological charge was interpreted as the baryon number. 

The model has been directly generalized to the case of the SU(3) group and subsequently 
to the generic case of the SU(N) groups [I]. Both SU(2) and SU(3) versions of the model have 
been canonically quantized using the collective coordinate formalism in [5J [6]. It was shown 
that the procedure of the canonical quantization leads to the appearance of new terms in the 
explicit form of the Lagrangian of the model that are interpreted as the quantum corrections to 
the mass of the skyrmion ('quantum mass corrections'). These quantum corrections restore the 
stability of the solitons that is lost in the semiclassical approach. Furthermore, this method of 
the canonical quantization has been subsequently generalized in [3 18] to the cases when the field 
U(x,t) belongs to a generic representation of the SU(2) and SU(3) groups and the stability of 
the solitons was explicitly shown. Interestingly, it appeared that the aforementioned quantum 
corrections are representation-dependent. 

The semiclassical quantization of the SU(3) Skyrme model has several shortcomings. For 
example it leads to a spectrum of masses of the baryon octet and decuplet and some physical 
characteristics of these that are not in close agreement with the values observed experimentally. 
One of the reasons for this disagreement is that the semiclassically-treated SU(3) model does 
not possess stable (semiclassical) solitons. Henceforth the classical solution of the SU(2) model 
(classical profile function) or some modification of it is used instead (see e.g. the overview [9]). 
An alternative way to overcome these problems is to consider the bound-state approach to the 
Skyrme model (see e.g. [TUl ITT] ). 

The purpose of the present paper is to show that the quantum mass corrections of the 
skyrmion that appear in the canonically quantized model are essential in ensuring the stability 
of the quantum solitons of the SU(3) model and realize Skyrme's original conjecture that 'the 
mass (of the meson) may arise as a self-consistent quantal effect' [2j. We find the stable quantum 
solitons by varying the complete quantum energy functional and then solving it numerically. 
The stability is ensured by the consequence of an iterative calculations. 
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Even though the SU(3) symmetry is not an exact flavour symmetry, by properly choosing 
the parameters of the model we obtain a baryon mass spectrum that is very close to the exper- 
imental one. We also focus on the influence of the Wess-Zumino-Witten (WZW), symmetry 
breaking (SB) and the quantum mass correction terms to the baryon masses and stability of 
the solitons. 

The paper is organized as follows. A brief description of the SU(3) model is given in the 
section 2 below. In section 3 we construct the quantum Skyrme model ab initio using the 
collective coordinate formalism. In section 4 the quantum energy functional is derived together 
with the asymptotic expression of its variation. Section 5 contains the numerical calculations of 
the mass spectrum of the baryons. Sections 6 accommodates the discussion and the concluding 
remarks. 

2 The setup of the model 

The Skyrme model is defined by the chirally symmetric Lagrangian density 

£ S k = -^Tr{R M R^} + ^Tv{[R fM K u ][K^R v ]} , (2.1) 

where the right chiral current is defined as R M = (d^U) J7L The pion decay constant f w and 
the dimensionless parameter e are the only parameters of the model. 

The main ingredient of the model is the unitary field U := U(x,t) that in addition to the 
fundamental representation (1,0) may also be defined for a generic irreducible representation 
(irrep) (A, /j.) of the SU(3) group. Then the basis states of the irrep (A, /jl) are labelled by the 
parameters (z,j,m) that are related to the hypercharge as y = |(// — A) — 2z (see [8] for the 
details). In such a way the classical SU(2) solitonic solution of the hedgehog type, that is 
defined by the canonical SU (2) SU (3) embedding, takes the following form, 

exp (i(a ■ x)F(r)) ^ U (x, F(r)) = exp (2i (j^ ■ x) F(r)) , (2.2) 

where F(r) is the soliton profile function, a are Pauli matrices and x is the unit vector. The 
generators J^\.\ represent the SU(2) subset of the SU(3) algebra generators J^imy 
superscript denotes that they are tensors of the adjoint representation (1,1) and thus can be 
expressed in terms of the Gell-Mann generators (we again refer the reader to [8] for the details) . 

As was shown in [8], the Lagrangian of the model depends on the irrep the unitary field 
U(x, t) was defined for. Interestingly the dependence on the irrep appears as an overall factor 
of the Lagrangian (|2.ip and is expressed in terms of the dimension of the chosen irrep and 
the eigenvalue of the quadratic Casimir operator of the SU(3) algebra. Likewise the baryon 
number of the model includes the same overall factor. Thus the Lagrangian of the model may 
be normalized in such way that at the classical level it is irrep-independent. However this is 
not the case at the quantum level. The canonical quantization of the model leads to a quantum 
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mass correction that is representation dependent. Furthermore, the Wess-Zumino-Witten and 
the symmetry breaking (mass) terms depend essentially on the chosen irrep [8]. In this work 
we shall restrict to the fundamental representation of the SU(3) group only. Nevertheless 
we shall be using the general formalism of [8], hence the generalization to the higher reps is 
straightforward. 

The hedgehog ansatz (|2.2p reduces the Lagrangian density (|2.ip to the following simple 
form, 

W» = = -{f (F'S + ^F) + ^^(2F'2 + (2.3) 

which represents the mass of the classical spherically symmetric soliton. Variation of this ex- 
pression leads to the differential equation for the profile function F(r) with topological boundary 
conditions F(0) = tt and F(oo) = 0. 

For completeness we also give the explicit expressions of the symmetry breaking and the 
WZW terms. The SU(3) chiral symmetry breaking term is defined as [8] 



-Csb = — Msb = 



2 r 

2 



Tr{[/ + [/t_ 2 .l}+2m2 Tr {([/ + [/+) J$$> Q) } ] . (2.4) 



The WZW action is given as an integral over the five dimensional manifold M 5 the boundary 
of which is the compactified spacetime, dM 5 = M 4 = S 5 x S . The standard form for this 
term is 

Swz(U) = J d 5 xe^Tv{K^K u K x K p K a } , (2.5) 

M 5 

where iV r is the number of colours. 



3 Canonical quantization 

The canonical quantization of the model is performed in terms of the collective coordinates. 
This approach allows the quantum unitary field to be cast in a factorisable form with the spatial 
and temporal dependent parts of the field being explicitly separated, 

U(x, F(r), q(t)) = D^(q(t)) U (x, F(r)) D^°\q(t)) . (3.6) 

Here D^'°\q(t)) is a Wigner D-matrix and is defined on the seven-dimensional homoee- 
neous space SU(3)/U(1), which is specified by the seven real, independent parameters q a (t), 
'the collective coordinates'. The ansatz (|3.6p may be effectively understood as the rotation of 
the field U$ in the quantum internal space parametrized by the collective coordinates q(i). 

The Lagrangian (|2.ip is considered quantum mechanically ab initio. Thus the collective 
coordinates q a and the conjugate momenta pp are required to satisfy the canonical commutation 
relations \pp,q a ] = —iS/3 a . On the other hand, this means that the coordinates q a {t) and the 
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velocities 4rq a (t) = q a (t) do not commute, rather they should satisfy the following commutation 



relation, 



-if a ^(q), 



(3.7) 



where f a ^(q) is a function of q a only and the explicit form of which will be determined by the 
consistency conditions below. 

The time derivative is defined by employing the usual Weyl ordering, 
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d 



doG(q) = -\q a , ^G(q) 



(3.8) 



where { , } represents the anticommutator. The ordering of the operators is fixed by the initial 
form of the Lagrangian (|2.ip . This allows us to avoid further ordering ambiguities in the case 
of the time derivatives (|3.8|) . The ansatz (|3.6|) is then substituted into the Lagrangian (j'2.1|) and 
followed by integration over the spatial coordinates. In such a way we obtain the Lagrangian 
cast in terms of the collective coordinates and velocities. Then the canonical momenta may 
be derived by restricting to the consideration of the terms of second order in velocities (this is 
because the terms of the first order in velocities vanish identically). Therefore the Lagrangian 
at the quantum level becomes 
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^Sk = r q a g a p{q, F) q p + [(q f - order terms] . 



Here g a p(q, F) is the metric tensor of the system and is expressed as 



9adq, F) = -C^ 7 - M )(g)(-l) z+M a J ( J P) * Z) _ Z , S I>P 5 M ^ M , C 



i(z',r,M') 



(<?), 



(3.9) 



(3.10) 



where (q) are functions of the coordinates q k only and the explicit form of which depends 
on the chosen parametrization of the SU(3) group. However the explicit form of Ca A \q) does 
not appear in the calculations. The quantum moments of inertia of the soliton are given by 
the integrals over the dimensionless variable r = ef n r, 



a,(F) 
ai(F) 



1 



1 



H (F) 
MF)- 



^—2tt [ dfr 2 (l -cosF) 
i J-K J 



1 



1 



F' 2 + 



1 

2P 



sin 2 F 



1 



dff 2 sin 2 F 



1 + F' 2 + 



1 



sin 2 F 



(3.11a) 
(3.11b) 



^ fw 3 

Note that ao(F) = and the summation in (|3.10p is over the basis states (Z,I,M) of irrep 
(1, 1) excluding the state (0, 0, 0). The quantum 'moment of inertia' a±(F) of the SU(3) model 
coincides with the quantum momentum a(F) of the SU(2) model. It is important to note 
that ai (F) is not equal to the mechanical momentum of inertia of the mass distributed by the 
classical spherically symmetric hedgehog field defined in (j2.3j) . 

The canonical momentum, which is conjugate to q@ is defined as 

dL S k 1 



(0) 



8qP 



{q a ,g a p} 



(3.12) 



4 



The superscript (°) was introduced to denote the canonical momentum obtained from (|3,9p . As 
we will show later, the WZW term shall contribute to the final form of the canonical momentum. 



Next, by requiring the canonical commutation relations 



(o) 
Pp > 



-id a R to be satisfied, the 



initially undetermined commutation relations fj3.7f) are constrained to be 



q ,<f 



ig af5 (q,F), (3.13) 



where g a ^(q, F) is the inverse of (|3,10p . This relation allows us to determine the explicit form of 
the (g)°-order terms in (|3.9|) . Thus after substituting (|3.6p into (|2.1|) . carefully manipulating the 
non- commutative variables and integrating over the spatial coordinates the additional quantum 
mass corrections are revealed [8]. Their explicit form will be presented in the section below. 

The contribution of the WZW term to the effective Lagrangian of the Skyrme model in 
the framework of the collective coordinate formalism was considered in [12] . By plugging (13.61) 
into (|2.5|) and employing Stokes's theorem and performing careful calculations, the WZW term 
takes the following form, 

Lwz(q,q) = -X' l -{q a ,C^(q)}, (3-14) 

where A' = ^= and B is the baryon number. 

The Lagrangian of the system with the inclusion of the WZW term becomes L' = Lgk+Lwz- 
The WZW term may be considered as an external potential of the system [13] . Therefore it 
shifts the canonical momenta pL°^ (|3.12|) by 



The metric tensor g a p and the functions f a " are not modified and the canonical commutation 
relations are preserved. 



4 The Hamiltonian 

The Lagrangian L' = ^sk+^wz effectively describes a system on a curved space with the metric 
g a p(Q>F) defined by (|3.1Up . The Hamiltonian for such a system is obtained by employing the 
general method of quantization on the curved space developed by Sugano et al. [T4]. This 
ensures the consistency of the Hamiltonian with the Euler-Lagrange equations of the model. 
We start by introducing seven right transformation generators 

R(A) = \{p« + XiCf\q), C[ a A) (q)} (4.16) 

that satisfy standard commutation relations of the SU(3) algebra. Here index A denotes the set 
(Z, /, M) excluding the case (0, 0, 0) , and C,%(q) are the reciprocal functions to Ca A \q), thus 
satisfy the standard orthogonality conditions. The generators -R(o,i,-) form a SU(2) subalgebra 
of SU(3) and may be interpreted as spin operators. This is because their action on the unitary 
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field can be realized as a spatial rotation of the skyrmion only. Next, it is convenient to define 
the eight transformation generator as -R(o,o,o) = or equally Ys = 1 in (|4.20p . In a similar 
way eight left transformation generators may be introduced, 

1 



L 



(B) 



R(A^,D 



t(1,1) (a)\ 



(4.17) 



using which the effective Hamiltonian of the model (with the constraint i?(o,o,o) = included) 
is found to be (see [8] for the details) 



H' 



1 t. ., Af f x/2 N 1 / 1 



2ai(F) 



(-l)AL {A) L ( _ A) -\"\+- 



2 I ai(F) ai(F) 



+ A Mi + AM 2 + AM 3 + M cl , 
where the following notation has been introduced: 



,1,—rn) 

(4.18) 



AMi 
AM 2 

AM, 



2vr 



a\{F) 



7T 



a\{F) 



r 2 dr sin 2 F 



sin 2 F 



(4.19a) 



r 2 dr (1 — cos F) 



,/r(2 - cos / ) + ^ ( ( 2 + cosF)jF '2 + 



ai(F)a 



2vr /" 



r 2 dr sin 2 F 



1 

2? 



/ 2 + A(f /2 + 



sin 2 F 



(4.19b) 
(4.19c) 



These negative quantum mass corrections appear because of the non-trivial commutation re- 
lations of the quantum coordinates and velocities (|3.13p . The kinetic part of the effective 
Hamiltonian is a differential operator constructed from the SU(3)/SU(2)-left and SU(2)~right 
transformation generators, thus the eigenstates of the model are 



\z T (Y T ),T,M T \ Z S (Y S ),5,M S ) ~~ V^ m (^ M ) D (Z t ,T,Mt)(Z s ,S,Ms)^ ^ ' 



(4.20) 



where the quantity D* on the right-hand side is the complex conjugate matrix element of the 
Wigner D-matrix for the (A, M) irrep of the SU(3) group and is expressed in terms of the 
quantum variables q . The topology of the eigenstates can be non-trivial and the quantum 
states contain an eighth 'unphysical' quantum variable q . 

Finally we are ready to consider the symmetry breaking term which takes the following 
form, 



L S B = -M SB = 4vr/ 2 J r 2 dr(l - cosF) 



1 



may 



(4.21) 



where the parameters m 2 , and m| are considered as the phenomenological parameters of 
the model. The expression (|4.2ip contains the operator Dv^j^Jq), which is a function of 
the quantum variables q a and acts non-diagonally on the states (|4.20j) . This means that 
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L^ z i Msb 7^ 0. Therefore the physical states of the system with the symmetry breaking 
term included need to be calculated by diagonalizing the total Hamiltonian [T5]. However the 
contribution of the symmetry breaking term is minor compared with the rest of the Hamiltonian 
and thus may be considered as a first order perturbation. The matrix elements of the symmetry 
breaking operator can be expressed in terms of two SU(3) Clebsch-Gordan coefficients, 



(AV) | n (l,l) ^| (A, M ) 
i T , z s , 

(\H) (1,1) (AV) 7 

Y t ,T,Mt Y^,T',M' T 



/ (AV) | D (1.1) (a \\ (A,m) \ _ 

\Y^,T',Mt, Ys' S ' > m 's' \o),(o) {q) \ Y T,T,M T ; Y s ,S,Ms> ~ 

r\ .A n -n f\> ,./^ 

(4.22) 



dim(A, (j,) ^ 
dim(A',/i') ^ 



(A,a«) (1,1) (A'y) 7 
y s ,5,M s Y^,S',M' S 



In the semiclassical approach the unitary field U(x,t) can be expanded in power series 
around the classical vacuum J7 = 1. In such expansion the parameters of the symmetry 
breaking term are obtained to be m\ = | ^m 2 + 2^m 2 ^\ and m| = ^= (^f 771 ^ ~~ m ^ where 
the experimental ratio ip- = 1.197 is imposed in order to obtain the standard mass terms of 
the 7T and K mesons. However we treat the model quantum mechanically ab initio and the 
collective coordinates q are not small perturbations. Thus the parameters m\ and m| need to 
be treated as generic parameters of the model. 

Putting all the ingredients together the energy functional of the quantum skyrmion in the 
operational form for the states in the irrep (A, M) becomes 
(^SU(3)/ i ]\/f'\ — A' 2 1 / 1 1 \ 

E(F) = 2 + 2 \ ajF) ~ ajF) J S < S + ^ + AM + Md + , (4.23) 

where AM = AM& and (Msb) represents the symmetry breaking operator Msb sandwiched 
between the states (|4.20p . The variation of the energy functional SE Jp"^ = gives an integro- 
differential equation for the profile function F(r) with the topological boundary conditions 
F(0) = 7r and F(oo) = imposed on top. At large distances this equation reduces to the 
asymptotic form 

f 2 F" + 2fF'-(2 + m 2 f 2 )F = 0, (4.24) 
where the dimensionless quantity fh? is defined as 

- 2 = - 4 (cf™ (A, M) - 8(8 + 1) - A' 2 + l) + 2 ^ + D + 3 

, 8AM 1+ 4AM 3 , AM 3 + 2AM 2 1 , , Mf x 

+ q- 7g\ + o~ + ~ (T?\~ fj?\ ) +\ M SB) ■ (4.25) 

3ai(Fj 2ai(F) ai(i<)ai(i 7 j ' 

2 2 

The tilded integrals AM/% and Msb are calculated using the dimensionless parameter r = ef n r. 
The quantity m = ef n rh is interpreted as the effective asymptotic mass of the baryon. For 
example, in the case of the nucleon it is = m n = 137.7 MeV. The corresponding asymptotic 
solution of (|4.24p is found to be 

/ ' m 2 1 \ 

F(r) = fe — + — exp(-mf). (4.26) 
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This solution is very important in ensuring the stability of the quantum soliton. It effectively 
translates into the requirement the integrals (|3,llal 13. 1 lb|) and AMk be convergent. Such a 
requirement is satisfied only if the asymptotic mass of the baryon m 2 > 0. This condition is 
only satisfied in the presence of the negative quantum mass corrections AM^ and symmetry 
breaking term (Msb) or at least one of them. However the general (non-asymptotic) integro- 
differential equation for the profile function F(f) obtained from the variation SE Jp = of 
the SU(3) model does not have stable solutions when the quantum mass corrections AM/% are 
absent. 

Finally we note that the symmetry breaking term is not necessary in ensuring the stability 
of the solitonic solution of the canonically quantized Skyrme model. The profile function F(r) 
has the required asymptotic exponential behavior (|4,26p even in the chiral limit when the 
symmetry breaking term is absent. In such way the canonically quantized Skyrme model is 
self-consistent. 

5 Numerical results 

We want to estimate the influence of the quantum mass correction AM on the stability of the 
quantum solitons and to compare the mass spectrum of the baryon octet and decuplet obtained 
using the semiclassical (rigid) and the quantum (soft) profile functions. 

Let us start by considering the semiclassical case first. The initial step is to find the classical 
profile function minimizing the energy functional of the classical SU(2) Skyrme model (|2.3p . The 
determined profile function asymptotically decays according to the power law F(f — > oo) ~ i 
and respects the topological boundary condition F(0) = it (see figure [1]). Then adding the 
symmetry breaking term modifies the profile function to be of the exponentially decaying form, 
F(f — > oo) ~ e~ m ^ r . The next step is to choose the parametrization scheme of the model. 
The SU(3) Skyrme model is parametrized by four parameters /„■, m 2 ,, rn\ and e. The first three 
parameters are of phenomenological origin, while the last one (e) is a dimensionless parameter 
that is usually constrained by requiring the model to fit the experimental data. Let us name 
these parameters the essential ones as they appear in the model explicitly. We shall also 
consider the following four phenomenological parameters: the nucleon mass = 939 MeV, 
the asymptotic nucleon mass m 3 ^ = m n = 137.7 MeV, the mean nucleon isoscalar (electric) 
radius (r 2 ) 1 / 2 = 0.78 fm and the mass of one of the heavier baryons (e.g. toa or ms) as possible 
input parameters for the model. We name them the fit parameters as they will be used to fit 
the model to the experimental data. 

The results of the numerical calculations using (|4.23p and based on the classical profile 
function are displayed in columns Mi and M2 of table [TJ Here the first column displays the 
experimental mass spectrum of the baryon octet and decuplet (the states are not discriminated 
by their spin polarization). The numbers standing at the right side of the mass show the devi- 
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ation (±%) of the calculated value from the experimental one. The parametrization of m\ and 



ml for column Mi is m 2 , = ± (m 2 + 2^m 2 K ) = 241 032 MeV 2 and mf = ^ \jrm 2 K 
384 638 MeV 2 , where = 137.7 MeV, m K = 495.7 MeV and / w = 92.2 MeV, 4£ = 1.197. 



MeVnxD 


Mi (%) 


M 2 (%) 


M 3 (%) 


M 4 (%) 


m as 


TV (939) 


939 (Input) 


939 (Input) 


939 (Input) 


939 (Input) 


137.7 (Input) 


A (1115.6) 


1015.0 (-9.0) 


1115.6 (Input) 


1167.0 (-4.4) 


1115.6 (Input) 


512.0 


S (1193.1) 


1091.0 (-8.6) 


1292.2 (8.3) 


1193.1 (Input) 


1202.1 (0.7) 


675.6 


S (1318) 


1129.1 (-14.3) 


1380.5 (4.7) 


1256.3 (-4.7) 


1235.9 (-6.3) 


741.0 


A (1232) 


1476.3 (19.8) 


1255.5 (1.9) 




1373.6 (11.5) 


377.7 


S* (1385) 


1523.8 (10.0) 


1365.9 (-1.4) 


1461.0 (5.5) 


1468.7 (6.0) 


496.3 


~* (1533.5) 


1571.4 (2.5) 


1476.3 (-3.7) 


1537.9 (0.3) 


1541.1 (0.5) 


590.2 


ft u (1672) 


1618.8 (-3.2) 


1586.6 (-5.1) 


1616.3 (-3.3) 


1600.7 (-4.3) 


670.2 


e 


5.7 


3.9 


4.5 


3.8 






Input 


60.9 


92.3 


58.3 




o 

m 


Input 


58 537 


(Input) 


211 224 




o 

m 8 


Input 


390 361 


546 807 


1 273 462 





Table 1: Baryon mass spectrum (MeV) obtained using (|4,23p . Columns Mi and M2 display 
the mass spectrum based on the classical SU(2) profile function. Column M3 is based on the 
quantum profile function of the nucleon. Column M4 is based on the individual quantum profile 
function for each state. The choice of the input parameters for each column is different and is 
emphasized in the table by (Input). The fourth input parameter for column M3 is the isoscalar 
nucleon radius (r 2 ) 1//2 = 0.78 fm. The third and fourth input parameters for columns M2 and 
M4 are the isoscalar nucleon radius and the asymptotic nucleon mass mf^ = = 137.7 MeV. 
The last column displays the asymptotic mass spectrum of the corresponding states obtained 
from the calculations of column M4. 

Note that the set of input parameters in both cases is different but always consists of both 
essential and fit parameters. The dimensionless parameter e is never an input parameter and 
is obtained by fitting the model to the experimental data. Let us explain both choices of the 
input parameters in detail. 

The standard choice in the semiclassical approach is to choose /„•, to 2 ,, mj and as the 
set of input parameters describing the model (column Mi of table [1]). However, restricting to 
the experimental value of f v even in the case of the SU(2) Skyrme model hardly reproduces 
the correct mass spectrum of the nucleon and its delta resonances [16]. Furthermore this 
choice leads to a value of the mean nucleon isoscalar radius, which may be evaluated using the 
following expression, 

(r 2 ) = 2 — I r 2 F' sin 2 Fdf , (5.27) 

ne z f4 J 

far from the experimental one. Thus there is no particular reason to restrict to this set of input 
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parameters and an alternative reasonable choice of the input parameters is m^, m^, m\ and 
(r 2 ) = 0.78 fm leading to a much better agreement with the experimental data (column M2 of 
table 

The approach we have been considering so far is not entirely semiclassical as we have been 
calculating the mass spectrum with the help of the classical profile function and (|4.23p which 
includes the quantum mass correction AM. However omitting this term leads to a complex 
value of the model parameter e and thus some other method to ensure the consistency of the 
model needs to be employed (see e.g. [9]). 

Let us turn now to the consideration of the self-consistent quantum SU(3) Skyrme model. 
The main difference with respect to the previous case is that instead of using the classical profile 
function we minimize the quantum energy functional (|2.3p by employing recursive calculations 
and thus obtain stable quantum profile functions for each state individually. The recursive 
calculations are performed in the following way: 

1. Find the classical profile function F m (f) minimizing the energy functional of the classical 
SU(2) Skyrme model (|2.3p and choose the set of input parameters describing the model 
as discussed above. 

2. Calculate the classical values of the integrals ai(F (0) ), ai(F (0) ) and AM(F (0) ) in (|3~2"3"j) 
and the (essential) model parameters by requiring the classical profile function to repro- 
duce the physical properties of the nucleon and arbitrary heavier baryon, e.g. A. 

3. Find the first approximation of the quantum profile functions F^\r) and '(f) by 
employing the asymptotic solution (|4.26|) and minimizing the quantum energy functional 
(|4.23p . i.e. solving the variational equation SE Jp"* = by using the classical values of 

the integrals ai(i ?(0) ), ai(-F (0) ) and AM(F m ), and the model parameters. Functions 

2 

Fff (f) and F^ (f) are found independently as they are describe the states with different 
quantum numbers. 

4. The obtained functions F$(r) and F^'(f) are used to calculate the updated values of the 
integrals a^Fff), ai(F$), AM(F$), and a^F^), ai(F^), AM(F^). The updated 
values of the model parameters are found by requiring the obtained profile functions to 
reproduce the physical properties of N and A. Then the procedure described in item [3] is 
repeated to get the second approximation of the quantum solutions Fj$ (f) and F^\f). 

5. The procedure described in item [J] is iterated until the convergent solutions Fn(t), 

F\(f) and stable values of the integrals ai(-Pjv), aif^jv), AM(Fn), and ai(F\), ai(F\), 

2 2 

AM (.Fa), and the model parameters are obtained. 

6. The obtained model parameters are used to find the quantum profile functions for the rest 
of the baryons. The same iteration procedure is employed (with the model parameters 
fixed) until the convergent solution and stable integrals are obtained. 
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F 




Figure 1: The classical profile function F(f) together with the quantum profile functions Fi{r) 
of the stable quantum solitons describing the baryon octet. The quantum profile functions 
correspond to the calculations presented in column M3 of table [TJ 

In case of the semiclassical approach this procedure fails - it does not lead to a stable 
soliton due to the absence of the quantum mass correction AM which not only contributes 
to the asymptotic mass of the state but also plays a crucial role solving variational equation 

SF ~ U - 

The mass spectrum of the quantum SU(3) model is presented in column M4 of table [TJ 
The choice of input parameters is the same as for column M2. Each state is described by an 
individual profile function obtained using steps 1-6 explained above and is displayed in figure 1. 
These quantum profile functions are very important as they can be used to calculate the 
magnetic moments and form factors of the corresponding states. The obtained mass spectrum 
is very close to the experimental one except for the A state. However delta resonances are not 
stable baryons; thus are not expected to be described by the model very well. 

The last column of table [TJ displays the asymptotic baryon mass spectrum. It reflects the 
mass density of the corresponding state in the asymptotic region, r — > 00. 

Finally let us discuss the calculations presented in column M3. The interesting fact is that 
such approach predicts a correct value of f n , while the approach used for columns M2 and M4 
leads to the value of f n much smaller than the experimental one. However this approach does 
not describe the A state as the corresponding integrals ()4.19p diverge. 
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6 Discussion 



In this work we have considered the stability of the topological solitons of the quantum SU(3) 
Skyrme model formulated in [8]. The model was shown to possess a family of stable quantum 
solitons whose energy functionals reproduce the mass spectrum of the baryon octet and decuplet 
in a good agreement with the experimental results. 

The semiclassical and quantum Skyrme models are essentially different models and lead 
to distinct integro-differential equations for the profile function F(f). In the semiclassical ap- 
proach the energy functional E(F) does not receive quantum corrections f|4. 19j) and the symme- 
try breaking term plays an important role in obtaining the exponentially decaying asymptotic 
profile function. Despite having correct asymptotic behavior, the semiclassical SU(3) Skyrme 
model does not support stable (quantum) solitons. Recursive solutions of the variational equa- 
tion SE J"p^ = do not converge and the classical profile function must be used instead. Hence 
the semiclassical Skyrme model is considered as describing a rigid quantum rotator because 
the profile function is fixed by the classical solution. 

The canonical quantization of the Skyrme model leads to the appearance of the quantum 
mass corrections ()4.19f) in its energy functional. These corrections not only ensure the correct 
asymptotic form of the profile function even in the absence of the symmetry breaking term, 
but also are necessary for obtaining stable quantum solitons. The recursive solutions of the 
variational equation SE Jp = with the quantum mass corrections present do converge and 
lead to the quantum profile functions which are different from the classical one. The difference 
is explicitly shown in the figure [1] where the classical profile function and the quantum profile 
functions for the baryon octet are displayed. 

Interestingly, the stability is preserved even if the Wess-Zumino-Witten and the symmetry 
breaking terms are not included in the model. Thus in this sense the quantum SU(3) Skyrme 
model is self-consistent and may be effectively understood as describing a soft quantum rota- 
tor. Despite the model being self-consistent the symmetry breaking term is necessary as it is 
responsible for the discrimination of the solutions with different hypercharges. Thus it must 
be included into the model in order to obtain physically reasonable results. 

The quantum approach to the SU(3) Skyrme model not only makes the quantum solitons 
stable, but also adjusts the model to fit better to the experimental results. Our numerical 
calculations of the mass spectrum of the octet and the decuplet of baryons presented in table [T] 
show that the quantum treatment of the model ah initio improves significantly the overlap 
with the experimentally-observed mass spectrum when compared with results obtained using 
the standard rigid rotator approach in the semiclassical version of the Skyrme model. The 
individual quantum profile functions obtained can be used to calculate the magnetic moments 
and form factors of the baryons. 
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